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Abstract 



We study the existence of a unique solution to semilinear fractional backward doubly stochastic 
differential equation driven by a Brownian motion and a fractional Brownian motion with Hurst 
parameter less than 1/2. Here the stochastic integral with respect to the fractional Brownian 
motion is the extended divergence operator and the one with respect to Brownian motion is 
Ito's backward integral. For this we use the technique developed by R. Buckdahn [3] to analyze 
stochastic differential equations on the Wiener space, which is based on the Girsanov theorem 
and the Malliavin calculus, and we reduce the backward doubly stochastic differential equation to 

r^. ' a backward stochastic differential equation driven by the Brownian motion. We also prove that 

the solution of semilinear fractional backward doubly stochastic differential equation defines the 
unique stochastic viscosity solution of a semilinear stochastic partial differential equation driven 

CN ' by a fractional Brownian motion. 
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1 Introduction 

This paper investigates semilinear fractional backward doubly stochastic differential equations (BDS- 
DEs) and semilinear stochastic partial differential equations (SPDEs) driven by fractional Brownian 
motion. Fractional Brownian motions (fBms) and backward stochastic differential equations (BSDEs) 
have been extensively studied in recent twenty years. However, up to now there are only few works 
that combine both topics. Bender [1] considered a class of linear fractional BSDEs and gave their 
explicit solutions. There are two major obstacles depending on the properties of fBm: Firstly, the 
fBm is not a semimartingale except for the case of Brownian motion (Hurst parameter H — 1/2), 
hence the classical Ito calculus which is based on semimartingales cannot be transposed directly to 
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the fractional case. Secondly, there is no martingale representation theorem with respect to the fBm. 
However, such a martingale representation property with respect to the Brownian motion is the main 
tool in BSDE theory. Hu and Peng's paper [9] overcame the second obstacle for the case of H > 1/2 
by using the quasi-conditional expectation and by studying nonlinear fractional BSDEs in a special 
case only. 

Nevertheless, there are many papers considering stochastic differential equations driven by frac- 
tional Brownian motion with Hurst parameter H > 1/2 ([2], [15] and references therein) or H < 1/2 
([13]), or covering both cases ([10]). For the case H < 1/2, one of the main difficulties is how 
to properly define the stochastic integral with respect to the fBm. In the paper of Cheridito and 
Nualart [6], and then generalized by Leon and Nualart [12], the authors have defined the extended 
divergence operator which can be applied to the fBm for H < 1/2 as a special case. In this paper 
we will use such definition for the stochastic integration with respect to the fBm, and then apply 
the non-anticipating Girsanov transformation developed by Buckdahn [3] to transform the semilinear 
fractional doubly backward stochastic differential equation driven by the Brownian motion W and the 
fractional Brownian motion B: 



rt rt 



Y t = i+ [ /(«, Y s , Z s )ds - f Z s I dW s + f j s Y s dB s , t e [0, T] 
Jo Jo Jo 



(1.1) 



into a pathwise (in the sense of fBm) BSDE 

Y t = i+jf( 3 Xes{Ts),Zse s (T s ))e- 1 (Ts)&s-fz s idW s , te[0,T]. (1.2) 

More precisely, the solutions (Y, Z) and \Y ,Z\ are linked together by the following relations: 
{(Y t ,Z t ),t€[0,T}}={(Y t {A t )e t ,Z t (A t )e t ^ ,t e [0,T]} 



and 



{(%Z t ) ,t€ [0,T]} = {(Y t {T t )et\T t ),Z t {T t )et\T t )),t€ [0,T]} 



where A t and T t are Girsanov transformations. It is worth noting that such kind of method was also 
used by Jien and Ma [10] to deal with fractional stochastic differential equations. 

It is well known that the solution of a BSDE can be regarded as a viscosity solution of an associated 
parabolic partial differential equation (PDE) (cf. [8] [18] and [23]), and the solution of BDSDE driven 
by two independent Brownian motions can be regarded as a stochastic viscosity solution of an SPDE 
(cf . [4] and [20] ) . So it is natural to consider the relationship between the solutions of our fractional 
BDSDE and the associated SPDE. We show that the solution of the above fractional BDSDE, which 
is a random field, is a stochastic viscosity solution of an SPDE driven by our fractional Brownian 
motion. To be more precise, the value function u(t,x) defined by the solution of a fractional BSDE 
over the time interval [0,£] instead of [0,T], see (4.4), will be shown to possess a continuous version 
and to be the stochastic viscosity solution of the following semilinear SPDE 

du(t, x) = [Cu(t, x) + f(t,x,u(t, x), V x u(t,x)a(x))]dt + 7tu{t,x)dB t , te[0,T], . . 

u(0,x) = $(x). [ ' 

Taking a Brownian motion instead of the fBm, equation (1.1) becomes a classical BDSDE, which was 
first studied by Pardoux and Peng [20]. The associated stochastic viscosity solution of SPDE (1.3) 
(with H = 1/2) was studied by Buckdahn and Ma [4]. Let us point out that, unlike [4] considering the 
stochastic integral with respect to B (H = 1/2) as the Stratonovich one and using a Doss-Sussman 
transformation as main tool, we have to do here with an extended divergence operator (H < 1/2), 
which condemns us to use the Girsanov transformation as main argument. However, this restricts 
us to semilinear equations. We will investigate the general case in a forthcoming paper, but with a 
different approach. 



The paper is organized as follows: In Section 2 we recall some preliminaries which will be used in 
what follows: Malliavin calculus for fractional Brownian motion, the definition of extended divergence 
operator and the Girsanov transformation. In Section 3 we prove existence and uniqueness results 
for stochastic differential equations driven by a fractional Brownian motion and backward doubly 
stochastic differential equations driven by a Brownian motion as well as a fractional Brownian motion. 
The relationship between the stochastic viscosity solution of the stochastic partial differential equation 
(1.3) driven by fractional Brownian motion and that of an associated pathwise partial differential 
equation is given in Section 4. 

2 Preliminaries 

The purpose of this section is to describe the framework that will be used in this paper. Namely, 
we introduce briefly the transformations on the Wiener space, appearing in the construction of the 
solution to our equations, some preliminaries of the Malliavin calculus for the fBm, and the left and 
right-sided fractional derivatives, which are needed to understand the definition of the extension of the 
divergence operator with respect to the fBm. Although the most of results discussed in this section 
are known, we prefer to provide a self-contained exposition for the convenience of the reader. 

2.1 Fractional calculus 

For a detailed account on the fractional calculus theory, we refer, for instance, to Samko et al. [25]. 

Let T > denote a positive time horizon, fixed throughout our paper, and let / : [0,T] — »■ K be 
an integrable function, and a € (0, 1). The right-sided fractional integral of / of order a is given by 

/?_(/)(*) = ^ [ j^-^ for a.a. , € [0,T]. 

Note that I^_ (/) is well-defined because the Fubini theorem implies that it is a function in L p ([0, T}), 
p>l, whenever / <G L p ([0,T]). 

We denote by I%_(L P ), p > 1, the family of all functions / e L p ([0,T}) such that 

f = Ix-(<P), (2-1) 

for some p E L p ([0,T}). Samko et al. [25] (Theorem 13.2) provide a characterization of the space 
I^_(L P ), p > 1. Namely, a measurable function / belongs to I^_(L P ) (i.e., it satisfies (2.1)) if and 
only if /eI"([fl,T]) and the integral 



s+e 



du (2.2) 



converges in L p ([0, T]) as e 4- 0. In this case a function ip satisfying (2.1) coincides with the right-sided 
fractional derivative 

(D T f)(s) = -i— (-J^r + a [ T f } a) ~ {M du) , (2.3) 

v T - JR ' r(l-a) \{T-s) a J s {u-s) 1 + a J 7 v ' 

where the integral is the L p ([0,T])-limit of (2.2). Moreover, it has also been shown in [25] (Lemma 
2.5) that there is at most one solution ip to the equation (2.1). Consequently, the inversion formulae 

I%_{DZ_f)=f, forall/el£_(27), 

and 

D T - (Ir-lf)) = f, for all / e L\[0,T\) 



hold. 

Similarly, the left-sided fractional integral and the derivative of / of order a, which are given, 
respectively, by 

W){x) = ^-.f o ^d, fora.a. *g[0,T], 

and 

( 2)« /)(a) = _i_fZW +a f Mlg dA (2.4) 

v °+- yA y r(l-a)V s Q Jo (s-w) 1+a / 

satisfy the inversion formulae 

/ « + (£>£+/) = /, for all / G I° + (27>), 

and 

^o+ ('?+(/)) = /, for all / G ^([0, T]). 

2.2 Fractional Brownian motion 

In this subsection we will recall some basic facts of the ffim. The reader can consult Mishura [15] and 
Nualart [16] and the references therein for a more complete presentation of this subject. 

Henceforth (17, J 7 , P) and W° — {W® : t G [0, T]} are the canonical Wiener space on the interval 
[0, T] and the canonical Wiener process, respectively. This means, in particular, that 17 — Co([0,T]) 
is the space of all continuous functions h : [0, T] i-» R with h(0) = 0, {W t °(w) = w(t), t G [0, T], cj G 0} 
is the coordinate process on 17, P is the Wiener measure on (17,23(0)) and T is the completion of 
£>(17) — o~{Wg,s e [0,T]} with respect to P. The noise under consideration is the process 



B t = I K H (t,s)dW°, t€[0,T}, 
Jo 



where Kh is the kernel of the fBm with parameter H G (0, 1/2). That is, 

H-l/2 r t 



K H (t,s) = C H 



(±\ ' [t - s) 11 - 1 ' 2 -{H- l/2)s l ' 2 - H f u H -V\u - s) H - l ' 2 &u 



where Cr = J (i-2H)8(i H 2H h+i/2) ■ The P rocess B — {B t : t G [0,T]} is an fBm with Hurst 
parameter H, defined on (17, J 7 , P) 1 i.e., B is a Gaussian process with zero mean and covariance 
function 

rsf\t -i 

R H (t,s):=E[B t B s ]= / K H {t,r)K H (s,r)dr = - (t 2H + s 2H - \t - s\ 2H ) , s,te[0,T}. 
Jo * 

Let T-L h be the Hilbert space defined as the completion of the space L e (0,T) of step functions on 
[0, T] with respect to the norm generated by the inner product 

(l [0 ,t],l [0 ,s]) nH = R H (t,s) = E[B t B s ], £,sG[0,T]. 

From Pipiras and Taqqu [24] (see also [16]), it follows that Hh coincides with the Hilbert space 

A^ /2 ~ ff := [/ : [0,T] -> R : 3 <p f G L 2 (0,T) such that /(«) = uVa-ff/V 2 -* (* ff-1/ V/(*)) («)} 

equipped with scalar product 

(f,g) A i /2 - H =C%V{H + l/2)*{<pf,< Pg ) L , ( p tT) . 



So l[o.t] i-» -B* can be extended to an isometry of A T onto a Gaussian closed subspace of 

L 2 (£l, T, P). This isometry is denoted by <p >-¥ B(ip). Moreover, by the transfer principle (see Nualart 
[16]), the map K : A^ /2 ~ ff -> L 2 ([0,T]), defined by (sec (2.3)) 

(10p)(8) = C H T(H + l/2)s l l 2 - H (py*- H u H -V 2 <p{u)) (s), a £ [0, T], 
is an isometry such that 

B(<p) = f (!C<p)(s)dW° and /Cl [0 , t] = K H (t, ■)l m , t € [0,T]. 



Using the properties of IC, Chcridito and Nualart [6] have extended the domain of the divergence 
operator with respect to the fBm B. This extension of the divergence operator in the sense of Malliavin 
calculus holds also true for some suitable Gaussian processes (see Leon and Nualart [12]). For B, this 
extension is introduced as follows. 

The following result identifies the adjoint of the operator K. (see [12]). It uses the left-sided 
fractional derivative Dq + defined in (2.4). 

Proposition 2.1. Let g : [0, t] — >• R be a function such that u i->- u 1 / 2 ~ H g(u) belongs to 



A/2-H 



(L«([0,a])) for some q> (1/2 - H)- 1 V H' 1 . Then, g £ Dom JC* , and for all u £ [0,T], 



{K*g){u) = C H T(H + 1/2)^-^ Dli 2 -" (s^^gis)) (u). 

Let iS (resp. <S/c) denote the class of smooth random variables of the form 

F = f(B(i Pl ),...,B( i p n )), (2.5) 

where ipi, . . . , ip n are in A T ' (resp. in the domain of the operator K*K) and / £ Cf (M n ). Here, 

Cf°(M") is the set of C°° functions / : R n — > R such that / and all its partial derivatives have 
polynomial growth. 

The derivative of the smooth random variable F given by (2.5) is the A T ' -valued random 
variable DF defined by 

DF = J2^(B(ip 1 ),...,B(< fin ))ip i . 

i—l % 

Now we can introduce the stochastic integral that we use in this paper. It is an extended divergence 
operator with respect to B. 

Definition 2.2. Let u £ L 2 (il, J 7 , P; L 2 ([0,T])). We say that u belongs to Dom 5 if there exists 
5(u) £ L 2 (Q) such that 

E [{K*KDF, u)z,2 ([0iT ])] - E [F5{u)] , for every F £ S K . (2.6) 

In this case, the random variable 8{u) is called the extended divergence of u. 

Remark 2.3. i) In [12] it is shown that the domain of JC*JC is a dense subset o/A T . Therefore, 
there is at most one square integrable random variable 6(u) such that (2.6) holds. 

(ii) In [6] and [12], it is proven that the domain of 5 is bigger than that of the classical divergence 
operator, which is defined by the chaos decomposition approach (see Nualart [17]). 

(Hi) In Section 3 and Section 4, we use the convention 



/ u s dB s = (J(ul [0)t ]), 
Jo 



whenever «l[o,t] G Dom S. 



2.3 Girsanov transformations 

In this section we introduce the Girsanov transformations on O, which we consider in this paper. 

In all which follows, we assume that 7 is a square-integrable function satisfying the following 
hypothesis: 

(HI) 7l [0it ] belongs to A^ /2 ~ H , for every t e [0,T]. 

We emphasize that Leon and San Martin [13] (Lemma 2.3) have shown the existence of square- 
integrable functions satisfying the above hypothesis. 

Now, for t € [0, T], we consider the transformations on Q of the form 



and 



T t (u)=u+ (/C 7 l [0;t] )(r)dr 



A t (uj)=u>- (/C7l[ ,t])(r)dr. 

Jo 



Notice that A t T t and T t A t are the identity operator of fl, and that the Girsanov theorem leads to 
write 

B(<p)(T t ) = B(<p) + f (£7l [ o,t])(r)(XV)(r)dr = B(p)+ / y r (K*K<p)(r)dr, 

Jo Jo 

for all ip £ Dom (IC*IC), and 
with 



E[F] = E[F{A t )e t ], 



e L = exp (J lr AB r -±J ((/C 7 l [0 , 4] )(r)) 2 dr") . 



We will need the following estimate of the above exponential of the integral with respect to the 
fractional Brownian motion: 

Lemma 2.4. Let 7 : [0,T] ^ M, 7 e L p [0,T] n Df[Q,T\, for some p > l/H, where D^[0,T] = {7 : 
[0,T] h-> M| J„ (J ip(x,t)dt) p dx < 00} and set ip(x,t) = !?_ vi/i-H ^{o<x<t<T}- Then there exists a 
constant C(H,p) only depending on H and p, such that 



E 



exp ^ sup 

0<0<T 



IsdB, 



< 2 exp 



|l/2(c(i/,p)G p (0,r,7) + 4V2) 2 |, (2.7) 



p\ 1 /p 



where G p (0,T, 7 ) := hll^o.r] -T H ~^ +T 1 '*- 1 /* (£ (/J ^(x, t)dt)") 

Proof: Let I7. = sup 0<t<T | J 7 s di? s |. According to Lifshitz[14] Theorem 1, P.141, and its corollary, 
for all r > 4\/2D(T, A/2), we have the inequality 



P{I't > r} < 2 1 - $ 



4 v / 2£>(T,A/2)* 



(2.8) 



where $(x) = -7= J_ exp{— y 2 /2}dy, D(T, A/2) is the Dudley integral (for more details, we refer to 



Mishura[15]), and A 2 = sup te r 0T i E 



IolsdB 6 



Since E [exp{/|,}] = 1 + J °° exp{x}P(Z£ > x)dx, by using the estimate of (2.8) we have 
E [exp{I£}] 

piV2D(T,\/2) /-H-oo 

=1 + / exp{x}P(Pj> > x)dx + / exp{x}P(Pj. > x)dx 

JO J±V2D(T,\/2) 

<exp{4V2^(T,A/2)} + 2 / + °° exp { x} f 1 - ft fizf^M^ d , 

*- J J4V2D(t,\/2) \ \ A yy 

= cxp J4y2L>(T, A/2) | + 2 / exp J4\^D(T, A/2) + xj [ 1 - -= / cxp{-y 2 /2}dy J dx 

<2cxp JA 2 /2 + 4V2D(T, A/2)| . 

Moreover, from Theorem 1.10.6 of Mishura[15] and its proof we know that 

\<C 1 (H,p)G p (0,T n ) 

and 

P(T,A/2)<C 2 (p)G p (0,T, 7 ). 

By substituting them to the former inequality, we easily get the wished result. ■ 

3 Semilinear fractional SDEs and fractional backward doubly 
SDEs 

3.1 Fractional anticipating semilinear equations 

In this subsection we discuss the existence and uniqueness of solutions to anticipating semilinear 
equations driven by a fractional Brownian motion B with Hurst parameter H € (0, 1/2). This type of 
equation was studied by Jien and Ma [10], and since it motivates the approach in our work, we give 
it in details. 

We consider the fractional anticipating equation 



X t =£+ [ b(s,X s )ds + [ ls X s dB s , te [0,T] 
Jo Jo 



(3.1) 



/() 

Here £ € L P (Q), p > 2, and b : Q x [0,T] xK^Risa measurable function such that: 
(H2) There exist v £ L 1 ([0,T]), v > 0, and a positive constant L such that 

\b{u,t,x)-b{w,t,y)\ < v t \x-y\, \ v t dt < L, 

Jo 

|6(«,t,0)| < L, 

for all x, y <G R and almost all w E fi. 
We observe that the above assumption guarantees that the pathwise equation 

( t (x) = x+ [ £ ; 1 (r ;s )6(r ;s , s ,£ s (r ;s )C s (x))d s , te [o,t], 

Jo 

has a unique solution. Henceforth, we denote it by £. 

Now we can state the existence of a unique solution equation (3.1): 



Theorem 3.1. Under Hypotheses (HI) and (H2), the process 

Xt = e t Ct{A t ,S(A t )) 



(3.2) 



is the unique solution in L 2 (Vt x [0,T]) of the equation (3.1) 7 such that 7Xl.r0 t i € Dom 8, for all 
t€[0,T}. 

Proof: We first show that the process X given in (3.2) is a solution of equation (3.1). For this we 
first observe that X belongs to L 2 (il x [0, T]) and we let F e Sic. Then, by the integration by parts 
formula and the Girsanov theorem, together with the fact that — ^-^- = j t (IC*ICDF)(Tt,t), we have 



E [FX t -F$=E [F(T t )Ct(0 - FCo(0] = E I — (F(T a )C.(0)d* 

= f ls E[(JC*JCDF)(T s , s )( s (0}ds+ f E[i?(T 8 )e7 1 (T0&(T a ,*,e 8 (T 8 )C.(0)]d« 

Jo Jo 

= I (i s E[(K,*1CDF)(s)X s }+E[Fb(s,X s )})ds. 

Jo 

Hence, since 7Xlr 0]4 i is square integrable, Definition 2.2 implies that 7-X"lrn t i belongs to Dom 8 and 
the equality in (3.1) holds, for all t e [0,T]. 

Now we deal with the uniqueness of equation (3.1). For this end, let Y be another solution of 
equation (3.1), F eS K undte [0,T]. Then, 



E[Y t F(A t )]=E[SF(A t )]+E 



F(A t )b(s,Y s )ds 



E 



ls Y s {K*KDF{A t )) (s)ds 



Therefore, the integration by parts formula, Fubini's theorem as well as the fact that — 4— — 
- ls {K,*KDF(A s ))(s) yield 

E [Y t F] 



=E[^F(A t )]-E 



C / j s (IC*ICDF(A s ))(s)ds 



E 
E 



F(A s )b(s,Y s )ds 



- E 



j s Y s {JC*JCDF(A s ))( s )ds 



-f r (!C*JCDF(A r ))(r) / b(s,Y s )dsdr 
I) Jo 

, * rr 
- E 



JO 



j r (JC*JCD{{IC*ICDF{A r ))(s)))(r)-f s Y s dsdr 



Hence, by using that Y is a solution of (3.1), Definition 2.2 and the relation 

{JC*JCD{{JC*JCDF{A r )){s))){r) = {K*K,D{{K*K,DF{A r )){r))){s), 
we obtain 



E[Y t F(A t )]=E[SF]+E 



F(A s )b(s,Y s )ds 



Consequently, by using the Girsanov theorem again, we get 



Y t (T t )ei\T t )=Z+ / 6(T s , S ,T s (T ;s )) es - 1 (T ;s )d S , 



which implies that Y t (T t )e t 1 (T t ) — C(£). That is, Y t — e t (t(A t ,t;(At)), and therefore the proof is 
complete. ■ 



3.2 Fractional backward doubly stochastic differential equations 

In this section we state some of the main results of this paper. Namely, the existence and uniqueness 
of backward doubly stochastic differential equations driven by both a fractional Brownian motion B 
and a standard Brownian motion W . 

Let {B t : < t < T} be a one-dimensional fractional Brownian motion with Hurst parameter 
H e (0,1/2), defined on the classical Wiener space (n',T B ,P B ) with ft' = C ([0,T];R), and 
{Wt = {Wf, W 2 , • • • , Wf) : < t < T} a d-dimcnsional canonical Brownian motion defined on the 
classical Wiener space (fi", F W ,P W ) with ft" = C o ([0,T]; R d ). We put (ft, J"°,.P) = (ft', T B , P B ) 
(Of 1 , F w , P w ) and let T = T° V TV, where TV is the class of the P-null sets. We denote again by B 
and W their canonical extension from ft' and ft", respectively, to 17. 

We let J^ = a{W T -W s ,t<s<T}vyV, J"/ 5 = a{B s ,0 < s < t} VVV, and & = J^VJf,ie 
[0, T]. Let us point out that .T 7 /^ is decreasing and F B is increasing in t, but <7t is neither decreasing 
nor increasing. We denote the family of er-fields {Gt}o<t<T by G. Moreover, we shall also introduce 
the backward nitrations H = {rl t = FYt v -^t }te[o,T] an d ^ W = {^tr}te[o,T]- 

Let 5^ denote the class of smooth random variables of the form 

F = f(B{ Vl ), . . . , B{<p n ), Wtyi), ..., W{1> m )), (3.3) 

where ipi, . . . ,ip n are elements of the domain of the operator K,*fC, ipi, . . . , ip m € C([0, T], R d ), / <G 
Cp°(R" +m ) and n,m > 1. Here, C£°(R" +m ) is the set of all C°° functions / : M" +m H- M such that 
/ and all its partial derivatives have polynomial growth. 

The Malliavin derivative of the smooth random variable F w.r.t. B is the Aj -valued randor 
variable D B F defined by 



^ dx, 

i—1 



D F = > . ^f( B (^l). ■ • ■ ' B( Vn ),W(^), ..., WW m ))<Pi, 



and the Malliavin derivative D F of the smooth random variable F w.r.t. W is given by 

Definition 3.2. (Skorohod integral w.r.t. B. Extension of Definition 2.2.) We say that u € L 2 (Q x 
[0,T]) belongs to Dom 8 if there exists a random variable 8 (u) E L 2 (ft) such that 

E [(IC*ICD B F,u) L 2 {[0tT]) ] = E [FS B {u)] , for all F € S' K . 

We call 8 B (u) the Skorohod integral with respect to B. 

Definition 3.3. (Skorohod integral w.r.t. W '.) We say that u e L 2 (ft x [0,T]) belongs to Dom 5 W if 
there exists a random variable 8 (u) G L 2 (fl) such that 



E 



f (DYF)u s ds 

Jo 



E [F6 w (uj\ , for all F E S' K 



We call 8 W (u) the Skorohod integral with respect to W . 

For the the Skorohod integral with respect to W we have, in particular, the following well known 
result: 

Proposition 3.4. Let u <G L 2 (Vt x [0,T]) be M-adapted. Then the ltd backward integral L u s J, dW 3 
coincides with the Skorohod integral with respect to W : 

u s ldW s = 8 w (u). 
o 



The extension to il of the operators T t and A t introduced in subsection 2.3 as acting over Q', is 
done in a canonical way: 

Tt(wW) := (Ttw'.w"). ^(w'V) : = (jW,w")> for (wW) efi = !J'x ft". 

We denote by Lg(0,T;R") (resp., L 2 I (0,T;R™)) the set of n-dimensional measurable random 
processes {<ft,t <G [0,T]} which satisfy: 

i) E [ j* T |pt| 2 dt] <+oo, 

ii) y> t is Qt- (resp., Ht~) measurable, for a.c. t € [0,T]. 

We also shall introduce a subspace of £ G (0, T; R n ), which stems its importance from its invariance 
with respect to a class of Girsanov transformations. Recalling the notation 



1^ '•= sup 

te[o,T] 



lsdB s 



from subsection 2.3, we define Lq*(0,T;R x R d ) to be the space of all G-adapted processes (Y,Z) 
which are such that 



E 



cxp{^} / (\Y t \ 2 + \Z t \ 2 )dt 
Jo 



< oo, for all p > 1. 



For the space L G '*(0,T;M x R d ) we have the following invariance property: 

Proposition 3.5. For all processes (Y, Z) G L G '*(0,T;M x R d ) we have: 
i) (Y u Z t ) ■= (YtiTJe^iTt^ZtiTJetHTt)) e L 2 G *(0,T:R x R d ) and 

it) (Y t ,Z t ) := (Y t {A t )s u Z t {A t )s t ) e L 2 G *(0,T;RxR d ). 

Proof: Since the proofs of i) and ii) are similar, we only prove i): 

For the case of (Y, Z), from the Girsanov transformation and Lemma 2.4 we have 



E 



expipl?} 



Y, 



dt 



E [cxp{^}(|Y t (T t )| 2 + \Z t (T t )\ 2 )ei 2 (T t )} dt 
E [cxp{pI^(A t )}(\Y t \ 2 + \Z t \ 2 )e- t Y \ dt 



<E 



cxp ^ p sup 

0<t<T 



IsdB, 



sup I {K, 1 l m ){s){K, 1 l [0A ){s)ds } f (\Y t \ 2 + \Z t \ 2 )e^dt 

0<t<T J Q J 



<CE 



cxp{( P + i)^} I (\Y t \ 2 + \z t \ 2 )dt 

Jo 



< + oo, for all p > 1. 
Hence, the proof is complete. ■ 

We now consider the following type of backward doubly stochastic differential equation driven by 
the Brownian motion W and the fractional Brownian motion B: 



Y t = £+ [ f(s, Y s ,Z s )ds - f Z s | dW s + f j s Y s dB s , t e [0, T] . 
Jo Jo Jo 

Here £ e L 2 (n, E^, P) and / : O" x [0, T] x R x R d ^ R is a measurable function such that 



(3.4) 
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(H3) i). f(-,t,y,z) is ^-measurable, for all t G [0,T], and for all (y,z) G R x R d ; 
(H3) ii). /(•,0,0)e J L 2 (fix [0,T]); 

(H3) iii). There exists a constant C G R + such that for all (y l7 Zi), (y 2 , ^2) € R X R' 
|/(i,j/i,zi) -/(*, 2/2,^)1 < C(|j/i - j/ 2 | + 1 21 -z 2 |), a.e.,a.s. 



Remark 3.6. Let us refer to some special cases of the above BDSDE: 

i) If 7 = 0, equation (3.4) becomes a classical BSDE [Pardoux and Peng [19]) with a unique 
solution (Y, Z) G L\ w (0, T; R x R d ); 



//£ G R <md / : [0,T] 



are deterministic, we can choose Z = and Y G 



L (Q! x [0,T]) wit/i 7Yl[o.t] G Dom{6 ), £ G [0, T], as £/ie unique solution of the fractional SDE 



Yt=£+ f(s,Y a ,0)ds + 







/ 7-^dB., iG [0,T], 

Jo 



which can be solved due to subsection 3.1. 

iii) Pardoux and Peng [20] considered backward doubly SDEs in the case that B and W are two 
independent Brownian motions and in a nonlinear framework. 

We let Cl' := {u/ G fl'\I^(uj') < 00}, which satisfies P B (fl') = 1 from Lemma 2.4. We first establish 
the following theorem: 

Theorem 3.7. For all ui' G Q! , the backward stochastic differential equation 

Yt(u/, •) = £+/ / (s,Y s {oj',-)s s (T s ,uj'),Z s (uj',-)s s {T s ,oj')) e^{T s ,uj')ds- f Z s (lu' ,■) I dW s , (3.5) 
te [0,T], has a unique solution (y(oj', •), Z(u', ■)) G C(0,T;lx R d ). 



Moreover, putting I Yt(u>' , •), Zt(uj' , •) ) := (0,0), for u>' G Q /c , the random variable lY t (u>',U>"), 

Z t (u',u"j\ is jointly measurable in (w',w"), a™d (?, ^) G L|*(0,T;M x R d ). 

Furthermore, there exists a positive constant C {only depending on the L?-norm of '£ and K/ylrg t i, 
L? -bound <?//(•, 0,0) and f/ie Lipschitz constant of f) such that, for all uj 1 G fi': 



£ 



IV 





2 /- T ^ 


2 


sup 


W 


•) 


+ 


2*w' 


•) 


d* 


te[o,T] 


Jo 





<Cexp{2/£(w')}. 



(3.6) 



Proofs We put F s (u/,y,*) - / (s, ye s (T s , uj 1 ), ze s (T 8 ,u')) eJ^T^u/), * € [0,T],(y,z) G R x R d , 

w' £ ff . Obviously, 

i) F s (-,y,z) is tj s -measurable and F s (u>', -,y,z) is .Fj^-measurable, ui' G f2'; 

ii) F s (ui' ,u>" ,y,z) is Lipschitz in (y,z), uniformly with respect to (s,ui',u)"); 

iii) \F s (co', -,0, 0)| <q/(.,a,0, 0)| exp{^(w')}. 

Using i 7 ^, equation (3.5) can be rewritten as follows: 

Y t (u') = £+ f F s (w',Y s (u',-),Z s (u',-)^ds- f Z B (u',-) I dW s , te[0,T], u/ G O'. (3.7) 

Step 1: We begin by proving the existence: From the conditions i)— iii) and standard BSDE 
arguments (see: Pardoux and Peng [19]) we know that, for all ui' G fl 1 , there is a unique solution 
I Y(ui', ■), Z(u>', ■) ) G L^ w (0, T; R x R d ). On the other hand, the joint measurability of F s with respect 
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to (w'.w") allows to show that, extended to Q'xQ" by putting (Y t (u',-),Z t (u',-)\ :=(0,0), u>' e fi' c , 

the process (Y,Z) is H-adapted. Let us show that (Y, Z\ € L^(0, T; R x R d ). For this end, it suffices 
to prove (3.6). 



Let a/ G O' be arbitrarily fixed. By applying Ito's formula to 



Y t 



Y, 



It follows that at u/, P^-a.s 



= 2Y t 



(F t (Y u Z t ) 



dt-Z t i AW t - 



we have at ui'. P^-a.s 



dt. 



Y t 



ds = r + 2/ Y s F s [Y s ,Z s jds-2 Y s Z s idW s 

ft 



<e 



Cx 



Y s 



2 1 



Z, 



+ C 2 \F S (0, 0) ¥ ) ds - 2 / F S Z S 1 AWs 



Hence, at u> , P -a.s., 



Y t 



2 1 ft 



z> 



ds<r 



Cx 



Y K 



C 3 e- s J - )d*-2 / Y s Z s idW s . 



Taking the expectation with respect to P , we notice that 



E v 



1/2' 



Y t (u',-)Z t (oj',-) 



dt 



< \E 



w 



sup 
te[o,r] 



Y t (u>',- 



1/2 



E 



»• 



1\ V2 



^(w',- 



dt 



< +oo. 



Consequently, F^ J Q Y s (oj', -)Z s (uj' , •) | dW s = 0, and by taking the conditional expectation with 
respect to T® , we obtain 



E 



w 



%(w',-) 



2 1 



Z«(w',-) 



d.s 



< # w [e] + / c^ 

Jo 



£(«',-) 



ds + C 4 exp{2/£(w')}. 



Thus, from Gronwall's inequality, we have 



P 



IV 



W,0 



< (P [C 2 ] + C 4 exp {27J (w')}) cxp{Cii}, t € [0, T] 



which, combined with the previous estimate, yields 



E 



w 



\f{ 


Yt(u',-) 


2 

+ 


Z t (w' 





) dt 


Jo V 










/ 



< Cexp{2Pf(w')}. 



Y t (w',-) 



by using 



In order to get the estimate (3.6), it suffices now to estimate E w sup te r 0T i 

equation (3.7), Burkholder-Davis-Gundy's inequality and the above estimate. 

To prove that (Y,Z) belongs even to Lg(0,T;R x R d ), we have to prove the uniqueness in 

IJ(0,I;Mx M d ). 
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Step 2: Wc suppose that I Y" 1 ,.^ 1 ) and (Y 2 7 Z 2 ) are two solutions of equation (3.7) belonging 

to L|(0,T;M x R d ) (Notice that L|(0,T;R X R d ) C I|,(0,T;R x R d )). Putting AY t = ^ - Y 2 and 
AZ t = Z\ - Z 2 , we have 





AF t 


Jo L 


*•( 


Y s \Zl)- 


f s (?': 


?2 


)'. 


As 


- f AZ s ldW s , 
Jo 


te 


[0,1] 




By applying Ito's formula to 


AY t 


2 

, we get that 


E 


AY t 


2 


Jo 


AZ S 


2 

ds 


= 2E 


/ A? s 
.Jo 


F s (Y s \Zj)-F s (Y 2 ,Z 2 y 


ds 














<E 


/' 

Jo 


2 


(Co + C 2 ) 


AY S 


2 1 

+ 2 


AZ S 


2" 


ds 


} 



and finally from Gronwall's lemma, we conclude that AY t = 0, AZ t — 0, a.s., a.e. 

Step 3: Let us now show that (9, z\ is not only in L^(0, T; R x M d ) but even in L|(0, T; R x R d ). 
For this we consider for an arbitrarily given r G [0,T] equation (3.7) over the time interval [0,r]: 



*7 = e 



Fjr/,zj)ds 



z;±dw a , t€[o,r]. 



(3.8) 



Let "H t r := J*^ V J 7 , 8 , t G [0,r]. Then H T = {"H[} te [ , r ] is a filtration with respect to which VF 
has the martingale representation property. Since Ft(y,z) is Qt- and hence also 'HJ'-measurable, dt 

a.e. on [0,r], it follows from the classical BSDE theory that BSDE (3.8) admits a solution (Y,Z) € 



L| T (0,r;lxf). Due to the first step this solution is unique in Lj| T (0,T;R x R d ). Hence, [Y u Z t 
(Yf,Zt\, dt a.e., for t < r. Consequently, (Y t ,Z t ) is %J"- measurable, dt a.e., for t < t. Therefore, 
letting t It we can deduce from the right continuity of the filtration ¥ B that ( Y, Z\ is G-adapted. 

It still remains to prove that (9, Z) G L^'*(0,T;M x R d ). 

Step 4: For the proof that [Y,Z\ belongs to L G '*(0,T;R x R d ) we notice that, by the above 



estimates and Lemma 2.4, 
E 



Y, 



Z, 



dt 



E 



cxp{pI^}E 



T( 


9 t 


2 

+ 


z t 


Jo V 









dt\T4 



exp{pl£} 

/o 

< J B[Ccxp{(2 + p)/y}] < oo, for all p>l. 

Hence, the proof is complete. ■ 

Corollary 3.8. The process lY,Z) given by Theorem 3.7 is the unique solution of equation (3.5) in 
4<*(0,T;RxM d ). 

Now we state the main result of this subsection: 

Theorem 3.9. 1) Let (9,2) G 4'*(0,T;lx R d ) be a solution of BSDE (3.5). Then 



{(Y u Z t ),t€ [0,T]} = {(y t (i4 t )e t ,Z t (^ t )e t ),te [0,T]} G 4'*(0,T;R x 
is a solution of equation (3.4) wit/i 7Ylr 0]t i G Dom 6 B , /or a^ £ G [0, T]. 
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2) Conversely, given an arbitrary solution (Y, Z) G L G '*(0,T;R x W i ) of the equation (3.4) with 
7^1[o.t] € -Dom 5 B , for all t G [0, T], the process 

{(%Zt) ,*e [0,T]} = {(r t (T t )er 1 m),^(T t )er 1 m)) ,* e [0,T]} G 4'*(0,T;Hx R d ) 
is a solution of BSDE (3.5). 

From the Theorems 3.7 and 3.9 we can immediately conclude the following 
Corollary 3.10. The solution of equation (3.4) in L G '*(0,T;R x M. d ) exists and is unique. 



Proof (of Theorem 3.9): We first prove that, given the solution lY,Z) of equation (3.5), (Y,Z) 

defined in the theorem solves (3.4). For this we notice that for F being an arbitrary but fixed element 
of Sjc, from Girsanov transformation and from the equation (3.5), it follows 



E [FY t -F£\=E 



F(T t )Y t - FY 



--E 



F{T t )Y + F(T t ) / F s (Y s ,Z s )ds- F(T t ) / Z s | AW S - FY 



We recall that E 
lt {K,*JCD B F) (t,T t ), we have 



F(T t ) f * Z s I dW s = E / * DY F(T t )Z s ds . Thus, from the fact that f t F(T ( 



E [FY t - F£] =E 



Y / ls (JC*>CD B F)(s 1 T s )ds 



+ E 
- E 



F(T S )F S (Y s , Z s ) ds+ / lr {K*lCD B F){r, T r )drF s ( Y s , Z s ) dt 



t r t r t 



D™F(T s )Z s ds+ / / D™j r (JC*!CD B F)(r,T r )drZ s ds 

JO Js 



Moreover, from Fubini's theorem, the definition of the Skorohod integral with respect to W, and from 
Proposition 3.4, we obtain that 



E 



JO Js 


t 

DYlr(IC*ICD B F)(r, T r )drZ s ds 


= f IrE 

Jo 


t 

IrE 


(IC*ICD B F)(r,T r ) f Z s idW s 
Jo 


dr. 



Df (K,*K,D B F){r, T r )Z s ds 



dr 



Thus, by applying the inverse Girsanov transformation as well as Fubini's theorem, we obtain 
E [FY t - F£] =E 

+ E , • . , ■ , ■ 

Jo V Jo v 

r t -I r ft 

E 



F{T s )F s (? Sl Z s )ds- f DYF{T s )Z s ds 

v ' Jo 

ls {JC*JCD B F){s,T s ) (% + f F r (? r , Z, ) (!;■ - / Z, | dlF, ) d« 



--E 



[ Ff(s,Y s ,Z s )ds -E [ DYFZ s ds +E f ls {rZ*rZD B F){s)Y s ds 
Jo J YJo J Uo 



where, for the latter expression, we have used that Y s is a solution of (3.5). Since Z G Lq* (0, T; M. x 

R d ) (c L G (0,T;Rx M. d )), it holds 



E 



D^FZAt 



= E 



F I Z s I dW s 

o 
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Consequently, E 



f Q {K*KD B F){s) ls Y 8 ds] = E [F {y t - £ - /„* /(*, Y„ Z s )ds + f* Z s | dW s } 



holds 



for all F e S' K . From Proposition 3.5 we know that both, [Y, Z\ and (Y, Z), belong to Lg (0,T;R x 

R d ). Consequently, F t - £- f* /(s, F s , Z s )ds + J * Z s 1 dVF s e L 2 (0, J", P), for all i e [0, T]. Moreover, 
using (3.6), 



E 



--E 



\lA[o,t]{r)Yr\ dr 



= E 



lA[o,t](r)Y r (A r 



e,„dr 



lA[o,t]{r)Y r e r (T r )dr 



< CE 



<CE 



T 2 

I 1 2 ^ 

|7 r l[o, t ](r)| sup Y r cxp{/y}dr 

ii re[0,T] 



<C / |7rl[0,t](0| dr < °°' 



lA[o,t](r)Y r exp{/£}dr 



<C / |7rl[o,t](r)| £[exp{CJ£}]dr 

Jo 



Thus, according to the definition of the Skorohod integral with respect to B we then conclude 7^1 [q^j G 
Dom 5 and 

f 7sY s dB s = Y t - U + I f(s, Y s ,Z s )ds - f Z s I dW s \ , t € [0, T}. 

Hence Y t = £ + /„* f(s, Y s , Z s )ds - J* Z s | dW s + f Q ls Y s dB s ,t e [0, T}. 

Let us prove now the second assertion of the Theorem. For this end we let (Y, Z) e L G '*(0, T;Ex 
M. d ) be an arbitrary solution of equation (3.4) and F be an arbitrary but fixed element of S' K . Then 
we have 



E[Y t F(A t )} 
=E[ZF(A t )]+E 
=h +I 2 -h + h 



F(A t )f(s,Y s ,Z s )ds 



E 



F{A t ) f 
Jo 



Z s I dW s 



E 



F(A t ) f j s Y s dB s 

Jo 



where, using the fact that 4j:F(A t ) = —^ t {JC*JCD B F(A t )){t) and Fubini's theorem, 



h =E [£F] - E 



£ / ~f s (JC*lCD B F{A s ))(s)ds 



h=E 



F(A s )f(s,Y s ,Z s )ds 



- E 



lr (JC*JCD B F(A r ))(r) / f(s,Y s ,Z s )dsdr 
o Jo 



From the duality between the ltd backward integral and the Malliavin derivative D w (recall that Z 
is square integrable), we get 



h=E 
=E 
=E 



f Z s DYF(A s )ds 
Jo 

f Z s DYF(A s )ds 
Jo 

f Z s D™F(A s )ds 
Jo 



-E 
-E 
-E 



Z S D™ 



lr {K,*K,D B F{A r ))(r)dr 

[ [ Z slr DY(IC*ICD B F(A r ))(r)dsdr 
Jo Jo 



j s (IC*ICD B F(A s ))(s) Z r idW r ds 

JO 



d.s 
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Moreover, from the duality between the integral w.r.t. B and D B (observe that 7FI r 0jt i £ Dom S l 



we obtain 



Ia =E 



--E 



ls Y 8 (K*KD B F(A s ))(s)ds 



j s Y s (IC*ICD B F(A s ))(s)ds 



-E 
-E 



~f r (IC*!CD B (!C*!CD B F(A r )(r))(s))y s Y s dsdr 



Jo 
t 



ls (K.*JCD B F(A s ))(s) I j r Y r dB r ds 



Consequently, using the fact that (Y, Z) is a solution of equation (3.4), i.e., 

Y t = £ + [ f(s, Y„ Z s )ds - f Z s 1 dW s + f ls Y s dB s , t £ [0, T], 
Jo Jo Jo 

we obtain that E[Y t F(A t )} = E[£F]+E [/„* F(A s )f(s,Y s , Z s )ds\ -E \jl Z 8 Df F(A s )ds 

by applying Girsanov transformation again and taking into account the arbitrariness of F £ <S^, it 



Therefore, 



follows that for all F £ SI 



/o 



E 



J Z s (T e )ej\T e )D™Fds =E F U - Y t (T t )s^ (T t ) + J /( S ,F s (T s ),Z s (r s )) £ ; 1 (T s )ds| 



Now, since according to Proposition 3.5, (Y t (Tt)e^ l (T t ), Z t (T t )e^ 1 (T t )) G 4'*(0,T;lx R d ), we have 
^(Ttjer 1 ^) -£-/„' /(a, n(T s ), Z s (T s )) £s - 1 (T s )ds G L 2 (ft, ft, P). Therefore, 

F t (T t ) £t - 1 (T t )=^+ / /( S ,K s (T s ),Z s (r s ))e7 1 (r s )d s - / Z 8 (T 8 )e7 1 (T 8 )4.dW„o.a., (3.9) 
Jo Jo 

forallie [0,T], which means that (Y, Z) := {y t (r t )e t _1 (r t ), Z^T^ 1 (T t ) , t £ [0,T]} is a solution of 
equation (3.5). Hence, the proof is complete. ■ 



4 The associated stochastic partial differential equations 

In this section we will use the following standard assumptions: 
(H4) i). The functions a : R d ->• R dxd , b : R d -> R d and $ : R d -> K arc Lipschitz. 



(H4) ii). The function / : [0,T] 



x r x i x 



1— >• R is continuous, /(i, -,-,•) is Lipschitz, uniformly 



with respect to t and /(•, 0,0,0) G L 2 (Q x [0,T]). 
We denote by (X^ x )a< s <t the unique solution of the following stochastic differential equation: 



dX\' x 

X$' x = x G R d . 



-b (Xl> x ) ds-a (X^ x ) idW s , s£ [0, t] , 



(4.1) 



We note that this equation looks like a backward stochastic differential equation, but due to the back- 
ward Ito integral, the SDE (4.1) is indeed a classical forward SDE. Under our standard assumptions 
on a and b, it has a unique strong solution which is {Jv^lsgjo.tj-adapted. The following result provides 
some standard estimates for the solution of equation (4.1) (cf. [21] Lemma 2.7). 

Lemma 4.1. Let X Lx = {X l s x , s £ [0, t]} be the solution of the SDE (4.1). Then 

(i) There exists a continuous version of X tx such that (s,x) h- > Xl' x is locally Holder (C a ' 2a , for 

alla£ (0,1/2)); 

(ii) For all q > 1, there exists C q > such that, for t, t' £ [0, T] and x, x' £ M. d , 



E 



sup \X 

0<s<t 



t,x\1 



<c q (i + \ x n, 



(4.2) 
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sup 

0<s<T 



A sAt _ A sAi' 



(4.3) 



< C q [(1 + \x\« + \x'\ q )\t - t'\t' 2 + \x- x'\ q . 

With the forward SDE we associate a BDSDE with driving coefficient /: 

Y t,x = $ ( X t,xj + r / ( r , X l r : x , Y?' x ,Z t T ' x ) dr - [ Z^ x 1 dW r + f lr Y^ x dB r , s e [0, t] . (4.4) 
Jo Jo Jo 

By Theorem 3.7 and Theorem 3.9, the above BDSDE has a unique solution (Y f ' x , Z t,x ) given by 

(Yt> x ,Zl> x ) = (Y t s ' x {A s )e Sl Z t f{A s )e s ) ,a£ [0,t], 
where, for all u/ € fi', P w -a.s., 

Y^ («', •) - * K' x ) + f / (r, *«•*, $*•* («', •) £r (T r , «'), Z^ x (a/, ■) e r (T r , u/)) e^(T r , u')dr 



Zl' x («', •) 1 dW r , s€[0,t}. 

(4.5) 

Pardoux and Peng [19] [20] have studied BSDEs in a Markovian context in which the driver 
F r (x,y,z) is deterministic; here, in our framework the driver 

F r (u', x, y, z) ■- f (r, x, ye r {T r , u'),ze r (T r , w')) e~\T rt w'), (r, x, y, z) e [0, T] x R d x R x M d , 

is random but it depends only on B and is independent of the driving Brownian motion W . In the 
following we shall define X*' x , yp and Z^ x for all (s,i) e [0,T] 2 by setting X*' x = x, Y*> x = F t *' x 
and Zf' x = 0, for t < s < T. 



We have the following standard estimates for the solution: 

Lemma 4.2. For allp > 1, there exists a constant C v € M + such that for all (t, x), (t' , x') <G [0, T] xM. d , 
J E Of, P w -a.s., 



E 



\v 



sup 

0<s<t 



Y^ {"',-) 



Z? («/,■) 



d.s 



p/2 



E 



w 



sup 

0<s<iA*' 



Y!<*{<J, ■)-¥!•* (a/,-) 



tAt' 



<C p (l + |s| p )exp{p/J(a/)}; (4.6) 



i\x' /, ,/ 



Zt' x (u>,-)-Zl> x (a/,.) 



d.s 



<C p exp{pI* (u/)} (|* - x? + (1 + |x| p + |x'n|£ - t'| p / 2 ) 



(4.7) 



Proof: Let us fix any u/ £ O'. For p > 2, by applying Ito's formula to | V^f ' ^ | ^ , it follows that, P w -&.s., 
for se [0,t], 



t,x\ \P 



v p(p-l) 



Y?' x (uj',-) P 2 ^' x («',-) 



l*(*r)l 



■P 



p-2_ 



F r ^ (a/, •) y,.*-* (a/, •) F r U, X^X'* («', ') , ^ (w', ■) dr + ^*' x («', ■) 1 dW r . 
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Let Q<s<t'<t<T. We take the conditional expectation with respect to Tlf t on both sides of the 



above equality, and we obtain 

p 



E 



--E v 



IV 



■w 



l*W)l p |J^ 



E 



■pE" 



IV 



pip - 1) 



Y^{u',-) P ~ 2 Zp*(uj',-) 2 dr\TW 



Y r t ' x (oj',-) P 2 ??>*(u/,-) 



xF r {oj , ,X t r ' x ,Yr t ' x ^',-),Z t r ' x (w , ,-))dr\J' t , t 



<E 



w 



t,x\\P I -rW 






-pE v 



Lii 



p-2 



2 



^P pV 



t,x|P 



y?' x (w',-) ^(o/,-) +C p e?(w') dr|J? 



-IV 

t',t 



Thus, from Gronwall's inequality and (H4) we have 



Y!>* {«/, •) 



<C P \E" 



our. I TT*^I P I TrVK 
Sup A r | \J- t 

0<r<s 



+ exp{pJ£(u/)} < C p (1 + |X S ^| P ) exp{p/£(u/)} 



Consequently, by using Doob's inequality, we get from the arbitrariness of p > 1: 

v 



E 



w 



sup 

0<s<t 



Y^ (u/, •) 



<C p (l + |af)exp{pJ£(u/)}. 



The first result follows from Burkholder-Davis-Gundy inequality applied to L Zl' x (u/, •) J, dW^ 

(see, e.g. [21]). 

Concerning the second assertion, without loss of generality, we can suppose t > t' . Let < s < 
t" < f. 

Using an argument similar to that developed by Pardoux and Peng [21], we see that, for some 
constants 6 > and C > 0, 



E w 

+ 9E W 
<CE W 
+ CE W 



y r <<* (u/, •) - y/-*' (J, ■) P '" 2 zr {J, ■) %>*' (J, ■) 2 dr|^ >t 



e-P(T r ,uj') 



\T W 
\ r t" ,t 



vt,x yt ,x 



t',2:' C, .' 



Y?>* (u/, •) - Y* ■* (u/, •) dr|^ 



r VV 






$ (*«•*) - $ (X* >* ) \f$ tt + Cexp{p/; («')} ( E 

V 



Y't.x -yt',x' 



2p 



<W, t 



1/2 



*?■* (a/, •) - 1? ■* (a/, •) dr|J? 



-w 



Consequently, from Gronwall's lemma and according to Lemma 4.1, 



and 
E w 



sup 

0<s<t' 



Y^ (c/, •) - *?■* (u/, •) 



F s <<* (a/, •) - Y^ (u/, •) 



<C n 



yt.x yt',x' 



exp {pi? (u/)} , < s < i' < i, 



< C p exp{p/£ (a/)} ((1 + |af + [z'n |t - t'| p/2 + \x - z'| p ) . 
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Finally, with the help of the Burkholder-Davis-Gundy inequality together with Lemma 2.4, we deduce 
that for all p > 2, there exists C p such that 



E 



w 



sup 

0<s<t' 



t' ,x' I , J 



y^ («', o - y: •* («', o 



p/2 



&,•*(<*/,■)-%'*'{«/,■) 



d.S 



<C p exp{pI^ (u/)} (Jx - a/| p + (1 + |af + |a/| p )|t - t'| p / 2 ) . 

The case p > 1 follows easily from the case p > 2. This completes the proof of the lemma. ■ 

We now introduce the random field: u{t,x) = Y*' x , (t,x) <G [0, T] x M d , which has the following 
regularity properties: 



Proposition 4.3. The random field u(t, x) is T? -measurable and we have Y*' x (u/ , ■)= V/' s (ui' , •) 
= u(uj', s, X*' x ), P w -a.s., < s < t < T, u>' e Q! . 

Proof: From Theorem 3.7 with terminal time t, we know that Y*' x is J 7 ]^ V J-^-measurable. Hence 
u(t, x) = Y t ' x is T^ t V J^-measurable. By applying Blumenthal zero-one law we deduce that u is 
J^-measurable and independent of W . The second assertion is a direct result from the uniqueness 
property of the solutions of (4.1) and (4.5) (cf. [8]). ■ 

Lemma 4.4. The process {Y*' x ; (s, t) E [0,T] 2 ,x £ R d } possesses a continuous version. Moreover, 
\u(t,x)\ < Cexp{/*}(1 + |x|), P-a.s. 

Proof: Recall that, for s € [0,t], 

Y*' x = $(X*' X ) + j S f (r,X t r < x ,?;' x er (T r ),2 t r ' x Sr (T r j) e-\T r )dr - J' Z*' x 1 dW r . 
Let < t < t' < T, x, x' e R d . Then by Proposition 4.3, we have, 

u(t,x)-u(t',x') = E 
and, thus, 



'y*' x -Y*'*'\T2 



\u(t,x)-u(t',x')\ P 
<CE 



E 



yt, X _ yZ , X | j.JU 



\/-t.x ^Crt' ,x' 



\n 



c 



E 



yt' .x' ^t'x' 



where, P-a.s., by Lemma 4.2, 

p 



E 



\/-t,x xrt ,x' 



\n 



yZ .X _ yZ .X | j.* 

< Ccxp{pl^} ((1 + \x\p + \x'\P)\t - t'\P' 2 



r'lP 



and 



< 



< 



Y, 



t' .x \}-t' ,x I ~r-B 

- *t> WT 



f \f (s,xt> x ' ,Yf x ' e s {T s ),zi' x ' £s{T s ))e:\T s )\&s\ 
J C(V(T s )/M,0,0)(l + 



\TS 



t'x' 



ds\T4 



<C(1 + \x'\ p )\t - t'| p exp{p/£} + C\t - t'\ p ' 2 E 
<C(l + \x'\P)ex V {pP T }\t-t'\Pl 2 . 



Y> 



Z 



t ,x 



As\H 



p/2 



19 



Consequently, for all (t,x), (t',x') <= [0,T] x R d , p > 1, P-a.s., 

|«(t,x) - u(t',x')\ < Cexp{/£} ((1 + |z| + \x'\)\t - if/ 2 + \x - a/|) . 

Hence, 

E [\u(t, x) - u(t', x')\ p ] < C p ((1 + \x\ + \x'\Y\t - t'y' 2 + \x- x'\p) , 

and Kolmogorov's continuity criterion gives the existence of a continuous version of u. ■ 

Henceforth we denote by C the second-order differential operator: 

£:=-tr(aa*(x)D 2 xx ) + b(x)V x , 
and we consider the following stochastic partial differential equations: 

du{t,x) = [Cu(t,x) + f (t,x,u{t,x)s t (T t ),V x u{t,x)o-(x)e t (T t )) s^{T t )] dt, te [0,T]; 

«(0, a;) = $(»). { ' 



and 



du(t, x) = [Cu(t, x) + / (t, x, u(t, x),V x ii(t, x)a(x))] dt + jtu(t, x)dB tl t G [0, T]; 
u(0,x) = $(x). 



(4.9) 



Our objective is to characterize u(t, x) = i^ ' x and u(t,x) — u{A tl t,x)£t — Y t ' x as the viscosity 
solutions of the above stochastic partial differential equations (4.8) and (4.9), respectively. 

Remark 4.5. In fact, equation (4.8) is a partial differential equation with random coefficients which 
can be solved pathwisely, and the equation (4.9) is a stochastic partial differential equation driven by 
the fractional Brownian motion B. 

First we give the definition of a pathwise viscosity solution of SPDE (4.8). 

Definition 4.6. A real valued continuous random field u : Q! x [0,T] x M n- M is called a pathwise 
viscosity solution of equation (4.8) if there exists a subset fl of Q' with P'(Q ) = 1, such that for all 
u>' G £1 , u(uj' ', •, •) is a viscosity solution for the PDE (4.8) at u>' ' . 

For the definition of the viscosity solution, which is a well-known concept by now, we refer to the 
User's Guide by Crandall et al. [7]. 

For the proof that u(t, x) is the pathwise viscosity solution of equation (4.8), we need the following 
two auxiliary results. 



Lemma 4.7. {Comparison result.) Let {HZ) hold. Let [Y 1 {u' , ■) , Z 1 {u' , -)\ and [Y 2 (u/ , •) , Z 2 {w' , ■) 

be the solutions of BSDE (3.5) with coefficients (£i,/i) and (^2,/2) ; respectively. Then, if £,\ < £2 and 
fi < f 2 , it holds that Y t l {uj', ■) < Y 2 {uj', ■), t G [0,T], P w -a.s., for all u/ G Q! . 

Proof: For the proof the reader is referred to the comparison for BSDEs by Peng [22], or also to 
Buckdahn and Ma [3]. ■ 



Lemma 4.8. {A priori estimate.) Let w' G Cl' and let (Y 1 ,^ 1 ) and \Y 2 ,Z 2 \ be the solutions of 

BSDE (3.5) with coefficients (£1 , f\) and {£2, f 2 ), respectively, and put 5Y (&/ , •) = Y 1 {uj' , •) — Y 2 [u>' , ■), 
S^ = e-e andSf s {uj\-) ^ (h-h) ^XH^^^s^s^^JUuj^^Ss^,^ e^iT^uj'). Moreover, 
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let C be the Lipschitz constant of f\. Then there exist /3 > C(2 + v 2 ) + /i 2 , v > 0, \i > 0, such that 
P w -a.s., for 0<s<T, and for all J G h' , 



E 



IV 



exp{/3(T-s)} SY a (u',-) 



ir 



E 



w 



sup 

0<r<T 



SY r (u',-) 



< E 



SZ r (uj' , 



exp{/3T} |^| 2 + ^ / cxp{/3(T-r)}|<5/ r ( W ',.)| 2 dr 
l 1 Jo 



dr 



<CE 



w 



m 2 + \6f r (u/,-)\'dr 



(4.10) 
(4.11) 



Proof: Let ui' G Cl' . By applying Ito's formula to 



SY s (u)',-) exp{/3(T — s)}, we obtain 



SY s (u',-) exp{/3(T- s )}+ / exp{/?(T - r)} [0 SY r (w', •) + 5Z r (u',-) 



dr 



|^| 2 exp{/3T} + 2 / r5y r (^',-)exp{/3(T-r)} /i [r,Y*(w', -)s r (T r ,uj'), Z\(J , -)e r (T r , u') 



+ 2 / exp{p(T-r)}8Y r (ui',-)8Z r (ui',-)ldW r 



- h [r^{w\-)e r {T r ^'),Z^{uj\-)e r {T r ,uj , ))]e- 1 {T r ,uj , )dr, 
and by taking the expectation with respect to P w on both sides we get, for v > 0, [i > 0, P-a.s., 



E 



IV 



exp{/3(T- s )} SY s (w',-) + / exp{/3(T-r)} [(3 SY r (u/,-) + 6Z r (u/,-) 



dr 



<E w [\5eeMPT}} 

exp{/3(T - r)} 



£ 



IV 



((2 + ^ 2 )c + m 2 ) a^V,-) 



c 



$Z r (wV) 



l*/W,-)| s 

M 2 



dr 



Finally, by choosing j3 > C(2 + v 2 ) + fi 2 , with i/ 2 > C, we obtain that P-a.s. 



E 



w 



cxp{/3(r- s )} £Y s (u/ 



< E 



\v 



xp{/3T}|^| 2 + / cxp{/3(T-r)W|5/,V,-)| 2 dr 



which is exactly (4.10). 

Estimate (4.11) can be proven by the arguments developed for (4.7). ■ 

Let us now turn to the solutions of our SPDEs. The next theorem is one of the main results of 
this section. For this let tt' — < uj' G f2'|u(u/, •, •) is continuous > and notice that in light of Lemma 

4.4, P B (n') = 1. 

Theorem 4.9. The random field u defined by u(ui',t,x) — Y t ' x (ui') for all uj' G O' is a pathwise 
viscosity solution of equation (4.8), where Y 1 ^^' ,-) is the solution of equation (4.5). Furthermore, 
this solution u(u)',t,x) is unique in the class of continuous stochastic fields u : fl' x [0,T] x M. d h- > M 
such that, for some random variable r\ G L a (T B ), 

\u(J, t,x)\< n{uj'){l + |x|), (t, x) G [0, T] x R d , P B {duj')-a.s. 

Remark 4.10. The uniqueness of the solution is to be understood as a P-almost sure one: let Ui(i = 
1,2) be such that u,(u/) is a viscosity solution of PDE (4.8) at uj' , for all u>' G fij. Then, by the 
uniqueness result of viscosity solution of deterministic PDEs (see: Pardoux [18], Theorem 6.14) we 
know that Ui(u/ , •) = U2(u' , ■), for all u>' G fl[ fl fl^- In particular, u\ — U2, P-a.s. 
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Proof of Theorem 4-9: We adapt the method used in the paper of El Karoui et al.[8] to show that u is 

a pathwise viscosity subsolution of equation (4.8). Recall that the set W :— {us' e Cl'\I^(ui') < +00} 

satisfies P B (£l') — 1. We work here on the set 0' := {u/ € Cl'\u(ui' , ■, ■) is continuous} which satisfies 

P B (0') = 1 in light of Lemma 4.4. 

Now, according to the definition of the viscosity solution, for an arbitrarily chosen w'ff!', we fix 

arbitrarily a point (t,x) £ [0, T] x R d and a test function ip e C^° such that (f(t,x) = u(u>',t, x) and 

<P > «( W V,')- 

For t <G [0, T] and ft, > 0, we have, thanks to the Proposition 4.3 and equation (4.5), 

u(J,t,x)=u{u J l ,t-h,X t ^ h ) + f F r (u,',X t r ' x ,Y?>*(w',-),Z t r ' x (u',-))dr- f £*•>', •) 1 dW r . 

We emphasize that for fixed u/, this BSDE can be viewed as a classical BSDE with respect to W and 
we recall that Y t ' x h (u/, ■) = u (to 1 , t — h, X t ' x h ) . Now for the fixed u/ £ il' , it holds that 

<p(t,x)<<p(t-h,X t t ? h )+ I F r L',X^,Y^{uj',-),Z^{oj',-))dr- f &*(<*/, •) I dW r . 



t-h 



t-h 



Let (Y '"' (u 1 , ■), Z ' a ' (u/, ■) J elg(i-/i,t;lx M d ) be the solution of the following equation eval- 
uated at to': for s £ [t — h, t], 



FfV,-) = V (t-h,X\\) + J' F r {u J ',X^,Y t ^ h (^,-),Z t ^ h (o J ',-)) dr 



-t.x.h 



(4.12) 



z; ' (u', •) i dw r . 



t-h 



-t.x.h , 



(4.13) 



From Lemma 4.7, it follows that Y t ' ' (u)', •) > f(t,x) — u{uj' ,t,x). Now we put 
Y s ^(uj',-) = Y t ^ h (uj',-)-ip{ S ,Xl' x )- [' G(J,r,x)dr 

J t-h 

and 

#>', •) = zf' V, •) - (V xW )( S , X?), 

where G(ui' , s, x) — d s <p(s, x) — Cip(s, x) —- F s (u)' , x,ip(s,x), V x tp(s, x)o~(x)) . From the equations (4.12), 
(4.13) and Ito's formula we have 



t-h 



F r (u',X t r > x , v (r,Xp x ) + Y r t > x (uj',-) + I G(uj',s,x)ds, 

t-h 



{V x ya){r, Xp x ) + Z^ x {w', ■) ) - {d rV > - C<p)(r, X l r x ) + G{lo' , r, x) 
Putting 



dr- / Z;' x (u>',-)ldW r 

t-h 



5{J, r, h) =F r I J, X**, <p(r, X^ x ) + J G(u>', s, x)ds, (V xV xr)(r, Xp x ) 

- (d rV - C^)(r,Xl- x ) + G{u' ,r,x), 

we have \5(u>', r,h)\ < k(u>') \X\: x — x\ , r £ [0, t], for some J 7 !? -measurable and P B -intcgrable n : il' 1— > 
R + . From the a priori estimate (4.10), it follows that 



E 



IV 



sup 

t-h<s<t 



WV) 



t-h 



znu',-) 



dr 



< CE 



tr 



\5(oj',r,h)\ dr 



t-h 



= hp(u>', h). 

(4.14) 
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where p(oj', h) tends to as h — >• 0. Consequently, it yields 

i-t 



E 



w 



t-h 



E w 



At 



h^fp(uj , h). 



(4.15) 



Furthermore, we have Yf' x {u') = E w Y{' x {oj' ', •) 
S'(u',r, h) = -{d r ip - &p)(r, X l r ' x ) + G(u/, r, x) 



fl h 6'{u/,r,h)dr 



where 



+ F r [co',X t r .' x , ( p(r,X t r ' x ) + Y r t ' x (uj',-) + / G(u',s,x)ds, (V x ^a)(r,X^ x ) + Z\: X (J , ■) 

Jt-h 



From the fact that / is Lipschitz and the estimates (4.14) and (4.15), we get 



#>') 



< E 



w 



t-h 



\6(«/,r,h)\+ *?'>', •) + Zp x (u',-) 



dr 



= hp(u> , h). 



Thus, from (4.13) (for s — t) and since Y t ' (u>', ■) > if(t, x), we obtain J G(u>' , r, x)dr > —hp(u>', h). 
Consequently, 4 J t _ h G(u>', r, x)dr > —p(ui', h). By letting h tend to 0, we finally get for ui' € 0', 

G(w', t, x) — df(p(t, x) — £ip(t, x) — F t (a/, x 1 ip(t, x), V x f{t, x)a(x)) > 0. 

Hence u(u>', t, x) is a pathwise viscosity subsolution of (4.8). The proof of u being a pathwise viscosity 
supersolution is similar. 

The proof of uniqueness becomes clear from Remark 4.10. ■ 

In analogy to the relation between the solutions (Y, Z) and (Y, Z) of the associated BDSDE and 
the BSDE, respectively, we shall expect that u(t,x) = Y^ x = Y^ x {A t )e u (t,x) € [0,T] x R d , is a 
solution of SPDE (4.9). This claim is confirmed by 

Proposition 4.11. Suppose that u,u are C 0,2 -stochastic fields over CI' x [0,T] x M. d such that there 
exist S > ond a constant C$ yX > (only depending on S and x) with: 



E 



\w(t,x)\ 2 +>+ / (\VMs,x)\ 2+5 + \D 2 x Ms,x)\ 2+5 )ds 



<C s , x ,te[0,T], forw = u,u. (4.16) 



Then u(t,x) is a classical pathwise solution of equation (4.8) if and only if u(t,x) is a classical 
solution of SPDE (4.9). 

Proof: We restrict ourselves to show u(t,x) solves equation (4.9) whenever u solves (4.8). For this, 
we proceed in analogy to the proof of Theorem 3.9. Let F be an arbitrary but fixed element of S'^. 
By using the Girsanov transformation we have 

E [u(t, x)F - u(0, x)F] = E [u(A t ,t, x)e t F - w(0, x)F] = E [F(T t )u(t, x) - Fu(0, x)} 

=E[F(T t )u(0,x) - Fu(0,x)} 
rt 



E 



F(T t ) f [£u(s, x) + f(s, x, u(s, x)e s (T s ), V x u(s, x)a(x)e s (T s ))e- 1 (T s )] ds 
Jo 



As in the proof of Theorem 3.9 we use the fact that ^F(T t ) = ~f t {lC*lCD B F){t,T t ) to deduce the 
following 



E[u(t,x)F-u(0,x)F] 



=E 



u(0,x) / ls {IC*ICD B F(T s ))(s)ds 
Jo 
t 



E 



7 r (JC*lCD B F(T r )){r)dr£u(s, x)ds 



+ E 
+ E 



F{T S ) [Cu{s, x) + /(*, x, u(s, x)e s (T s ),V x u(s, x)a(x)e s (T s ))sj\T s )] ds 

lr (JC*JCD B F{T r )){r)drf(s,x,u{s,x)e s {T s ) 7 \7 x u{ s ,x)a(x)e s {T s ))e- 1 (T s )di 



Js 
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Thanks to the assumption that u(t, x) is a pathwise classical solution of equation (4.8), we obtain 
E[u(t,x)F-u(0,x)F] 



=E 



F(T S ) [£u(s, x) + f(s, x, u{s, x)e s {T s ),\7 x u{s, x)a(x)e s {T s )) £ - 1 (T s )] ds 



E 



--E 



7s (/C ICD B F{T s )){s)u{s,x)ds 



F [Cu(s, x) + f(s, x, u(s, x), W x u(s, x)a(x))\ ds 



E 



j s (K.*1CD B F)(s)u{s,x)ds 



Consequently, 
E 
=E 



j s (lC*lCD B F)(s)u(s,x)ds 



F(„( i ,x)-»(o,x)-/[r„(„) + /(«,„(«),v I „(« W x))]d,) 



From the integrability condition (4.16) we know that 

u(t, x) — u(0, x) — [£u(s, x) + /(s, x, u(s, x), V x u(s, x)c(x))] ds e L 2 (il, F 1 P). 



Moreover, 7lr ,t]W <E L 2 ([0,T] x fl). Indeed, 



E 



/ \lA[Q,t]{r)u{r,x)\ 2 dr 
Jo 



\ lr l m (r)\ 2 E[\u(r,x)\ 2 }dr < C 5 , x [ |>l[cu] (r)| 2 dr < oo. 
o Jo 



By Definition 3.2 we get 



E 



F j s u(s,x)dB s =E Flu(t,x) — u(0,x) — / [Cu(s,x) + f(s,x, u(s,x), V x u(s,x)a(x))] ds) 



It then follows from the arbitrariness of F € S'^ that 

u(t, x) = &(x) + / [Cu(s, x) + /(s, x, u(s, x), V x u(s, x)a(x))] ds + / ■y s ii(s,x)dB s . 
Jo ~ Jo 

The proof is complete now. ■ 

Remark 4.12. The regularity ofu in the above proposition is difficult to get under not too restrictive 
assumptions (like coefficients of class Cf b , linearity of f in z). 

Remark 4.13. Notice that generally speaking, a continuous random field after Girsanov transforma- 
tion At is not necessarily continuous in t any more. We give a simple counterexample: 
Let < s < T be fixed and 



u(t, x) — 



(t - s)B t 



< t < s; 



(t - s)sgn(B s ), s<t<T. 



It is obvious that u(t,x) is F B -measurable and continuous in t. But after Girsanov transformation 
A t , it becomes 

u(t,x) = u(A t ,t,x)e t = { (' - S > ( B < - JoW [M )Mav) e u < t < s; 

{ (t- S )sgn(B s - f s (/C 7 % t ]Xr)dr)e t , s < t < T, 
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which is not continuous in t on 



J : inf I B, 

t£[s,T] 



(X/ylp,*]) (r)dr ] < < sup (B s - / (/C 7 l [0 , t] ) (r)dr ) 
/ te[s,T] V Jo / 



However, as we state below, the random field u has a continuous version in our case. To this end 
we need the following technical result: 

Lemma 4.14. Let 7 be such that (HI) holds. Then there exist positive constants C and q such that 
for all r, v, 1/ G [0, T], v < v' , 



(£7l [v',v] ) ( s ) (^7l [0,r] ) ( s )ds 



<C\v-v'\ q . 



Proof: We have 







(/C7l[„/ iV ])(s)(/C7l[o, r ])(a)ds 



2/ 



1/2 



< I ^ (/C 7 l Ku] ) 2 (s)ds 



1/2 



(^7![o,r]) 2 (s)ds 



<C / (IC 1 l [v , M y(s)ds 



where the last inequality follows from [13] (Lemma 2.3). Also, by the proof of Lemma 2.3 in [13] and 
using the notation a = 1/2 — H , we have 



' r- a y r 



as I r 7 r 



(/C7I [„,,„] )(s) 

=1 i">] (s) (^ (s) + f(f^o /, F~^ dr J " Ws) rxT^) X' (r'- 8 )i+° dr (4 ' 17) 

=/i(s) + / 2 (s). 
Now, applying [13] (Lemma 2.3) again, we obtain 

r-T 



h(s) 2 ds= l[v'.v](s) 1[0,«](«) 

Jo V 



7 (s) 



as" 



r "7 r 



T(l - a) y v (r - s) 



- sU+a 



dr 



d.s 



(4.18) 



<| V _^|(P'-2W|| ||2 < C | v _ v /|( P '-2)/ P ' 



On the other hand, for m = 1 + rj and g m = m/rj, with 77 small enough, we can write 



v —a 

r a y r 



(r-s)!+ 
<|w-i/| 1/<z ' 



-dr 



v r- ma \ lr \ m 
„/ (r - s ) m ( 1+a ) 



dr 



l/m 



-^'l 1 /^ 



<C|v-t/| 1/<? " 



, ( r — s \m(l+a) 



-dr 



l/r, 



(0)\m e -r, / (j - s )- m ^+ a ) r - ma+ ^{9 -r^^^drdO 



l/m 



<c| v _ t /| 1/9m (v'-s)- 27 ' /m 



(0)\™$-v / (r-s)- mQ+I '- 1 r- ma+, '(6i-r) ma - r '- 1 drd6' 



l/r, 
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Hence, [13] (Lemma 2.2) gives 



v —a 

r Q 7 r 



-dr 



/„, (r- S )i+« 

<q v - v / l 1/9m (v / -s)~ 27 ' /r ' 

<qv-v / | 1/< ' m (t; / -*)~ a ~ , ' /m 



6 7 (0)| m 0-V ma+n ( v ' - s)- ma+r '{9 - s)- l {9 - v') ma ~' n M 



l/m 



i^wr^v ma+, '(0 - s )- i+ma -vd9 



1/r. 



which, together with (4.17), implies, for p < 1/a, p' = 271 



(l-(ma-tj)p) 



and 1/p' + l/q' = 1, 



/ 2 (s) 2 ds 



<C\v-v'\^\ {v 1 - sf q ' { - a - r)/m) As 



|0 7 (0)| m 0-"(0 - s)- 1+mQ -"d0 



2p' 



d.s 



<C|u-t/| 2/ «™. 
Thus, we get the wished result. ■ 

Lemma 4.15. The random field u(t,x) := u(A t ,t,x)et, (t,x) G [0, T] x R d /ias a continuous version. 
Proof: In the following, for simplicity of notations, we put 



For < u' < v < T, we notice that (Y t ' x (A v ) ) is the solution of the BSDE 



Y^{A V ) = $ (x*' x ) + / / (9*'*'"'") e-^T^Jdr - / 2^(4,) | d^ r , (4.19) 

Jo Jo 

while (^(A,/)) is the solution of the BSDE 

Y^{A V ,) = $ (**•*) + f / (e^^>') e-^T.^Odr - f Z^(A V ,) | dW r . 

We set J^ = cxpj J Q (/C7l[ Oi „])(s)(/C7l[ 0;r .])(s)ds k and we observe that e^ 1 (r r A t ,) = e^T 1 ^),/". 
Moreover, equation (4.19) can be written as follows: 

Y S ^(A V ) =$ (X*' x ) + f/ (r,X«> x , (f r *^(^),Z^(^)) £r (T r )(J r ")- 1 ) e-\T r )r r dr 



Z^ x (A v )idW r , s£[0,t], 



2(> 



and a comparison with (4.5) suggests the similarity of arguments which can be applied. So, by a 
standard BSDE estimate, see, for instance, the proof of Lemma 4.2, we have 



E 



sup 

,0<s<T 



Y^(A V )-Y^(A V> ) 



<CE 



<CE 



dr 



f(@ t f^ v )e- 1 (T r A v ) - /(e*'*.»'» )s-\T r A v ,) 

\f(Q^ v ) (e-\T r A v ) - e-\T r A v ,)) \ p + | (/(0*™») - f(&™"'j) e^{T r A v ,) 



Recalling that e r 1 (T r A v ) — e r 1 (T r )J^ and applying Lemma 4.14, we get 



dr 
(4.20) 



\s-\T r A v ) - e-\T r A v ,)\ p = £ - p (T r ) 



TV TV 

tj r rJ r 



<Ce~ p (T r ) 



--Ce- p (T r ) 



(K/yl[o,t;])(s)(/C7l[ 0iT .])(s)ds- / (Kyl [oy] )(s)(Kyl [0tr] )(s)ds 



(4.21) 



(Kyi [ v >, v ] ) (s) (Kyi [0 , r ] ) (s)ds 



< Ce~ p (T r ) \v - v'\ 



r\pq 



On the other hand, 

f(Qt,X,V,V) _ f(Qt,X,V,v' 
<C 



e- p (T r A v ,) 



Y^ X (A V ) P + Zl: x (A v ) P ) \e r (T r A v ) - e r (T r A v ,)\ p e~ p (T r A v ,) 



=C 
<C 



Y^(A V ) + Z^(A V ) 



cxp ■ 



(Kyl [v , M )(s)(Kyl { Q yr] )(s)ds \ - 1 



(4.22) 



Y^(A V ) P + ^(A^^lv-vT. 



Plugging estimates (4.21) and (4.22) into the equation (4.20), Lemma 4.2 yields that 



E 



sup 

0<s<T 



<CE 



Y^(A V )-Y^(A V ,) 



f(@ t f^ v )(e- 1 (T r ))\ p + Y^(A V ) + Z^(A V ) 



dr 



<C(l + \x\ p )\v-v'\ 



ripq 



For the latter inequality we have used the following estimate of Z, which proof will be postponed until 
the end of the current proof. 



Lemma 4.16. There exists a constant C p such that E 



su Po<s<t 



< C p (l + \x\p). 



We continue our proof of Lemma 4.15: According to the proof of Lemma 4.4 we have 



E 



sup 

0<s<T 



< E 



y^ x (a v )-y^(a v ) 

1\ 1/2 



E 



E 



sup 

se[o,r] 



Y. 



t ,X X/"t,X 



U : 



O A,-. 



sup 

a£[0,T] 



yt ,x yt,x 



2„ 



(E [e- 2 (T v )]) 1/2 < C ((1 + |s| p + \x'\ p )\t t'\ p / 2 + \x x'\ p ) 
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Hence, by combining the above estimates we obtain 
e\ sup \y^ x \A v ,)^Y^ x {A v ) 

</"» [ JTJ vt' ,x ( /[ \ \}t' ,x' 



sup 

0<s<T 



sup 

,0<s<T 



rt ,X 



Y: •* (A v ) - Y^{A V ) 



Y^ X {A V ,)-Y^ X (A V ) + E 

<C ((1 + \x\ p + \x'\ p )(\t - t'\P/ 2 + \v- v'\ pq ) + \x- a/| p ) • 

Consequently, from the Kolmogorov continuity criterion we know the process {Yp x (A v ); s,t, v € 
[0, T],x £ M. d } has an a.s. continuous version. From Lemma 4.14 we have that e t is continuous 
in t. It then follows that u{t, x) — Y t ' x (A t )st has a version which is jointly continuous in t and x. ■ 

Proof of Lemma J^.16: For simplicity we suppose all the functions $, /, er, b are smooth. For the proof 
of the general case, the Lipschitz functions have to be approximated by smooth functions with the 
same Lipschitz constants. We define (Y_ t,x , Z_' x ) to be the solution of the equation 

yt,* = & { xl> x )V x Xl> x + [' [f' x {%)e-\T r )W x Xl' x + f y (K)Y!r X + f z (K)Z!r X ] dr - f £ x | dW r , 

Jo Jo 

(4.23) 

r t,x \rt.x 



(4.24) 



where S; = [r,X^ x , Y^ x e r (T r ), Z^ x e r {T r )\ and 



V x Xl' x = 1-1 V x Xi-*l/ (Xl< x ) ds - [ V X X<'V (X**) I dW s , r E [0,*]. 

J t J r 



With the arguments used in the proof of Lemma 4.2, we get 



E 



sup \Y_ 

0<s<T 



t.x\P 



T \ P/2' 

|Zt' x | 2 dsl 







<C p (l + \x\ p ). 



(4.25) 



By arguments which by now are standard, it can be seen that the processes X l ' x ,Y X ' X and Z f ' x are 
Malliavin diffcrcntiablc with respect to W, and thus, for s < u < t < T, 



D™YZ*=V{XF)Dyxl*. , „,, .,, 



D™Zl' x ldW r 



(4.26) 



dr. 



fx{^r) £ r ( T r)D u X r ' x + f y (E r )D u Y r ' x + f z (^ r )D u Z/ 
On the other hand, from 

Y t.x = yt,x + f f (r, X$' x , Y*' x e r {T r ), Z^e r (T r )) e~ l (T r )dr - f Z^ | dW r , 9 < s < u < T, 
we get 

fU 

D^Zl x idW r 



D^Yl- x = - Ztf 



fAKYr(T r )D™xr + fj*)DW + f'AKW&s 



(4.27) 



dr. 



From the above three equations (4.23), (4.26) and (4.27) and the relation 

D™(Xl' x ) = V x X t 8 ' x (V x X$' x )- 1 <r(Xf.> x ) 



(4.28) 
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(see: Pardoux and Peng [20]), we obtain that, for s < u < t < T, 

F*/(V X X^)-V(X^), 



LJ U L s 



Finally, from standard estimates for (4.24) and from the estimate (4.25), we get 



E 



sup 

0<s<T 



<E 



sup \Y_ t f\ p sup \V x X* t ' x \ P sup \(j{X t f)\ p 

0<s<T ' 0<s<T 0<s<T 



sup \Y_ 
0<s<T 



t,x\P<H 



l/9l 



E 



SUp V^X 
0<s<T 



t,x\-P92 



1/92 



E 



sup \a(X 

0<s<T 



t.X\\Pl3 



1/93 



<E 

<C p (l + \x\P), 
where l/o/i + I/92 + 1/<Z3 = 1, with qi,q2,qs > 1- The proof is complete. ■ 

The above proposition motivates the following definition: 

Definition 4.17. A continuous random field u : [0,T] x M. d x (!' h> R is a (stochastic) viscosity 
solution of equation (4.9) if and only ifu(t,x) = u(T t ,t 7 x)e^ (T t ), (i,x) <E [0,T] x R d is a pathwise 
viscosity solution of equation (4.8). 

As a consequence of our preceding discussion, we can formulate the following statement: 

Theorem 4.18. The continuous stochastic field u(t,x) := u(A tl t,x)e t = Y t ' x (A t )e t — Y t ' x is a 
stochastic viscosity solution of the semilinear SPDE (4.9). This solution is unique inside the class C.f 
of continuous stochastic field u : Q! x [0, T] x M n- M such that, 

\u{t,x)\ <Ccxp{/£}(l + |x|),(i,x) e [0,T] xM d ,P-a.s., 

/or some constant C only depending on u. 

Remark 4.19. 1) It can be easily checked that u e C p if and only if (u(At, t, x)et) € Cp if and only 

zf(n(T t ,t,x)ei 1 (T t ))eC p 3 . 

2) As a consequence of the preceding theorem we have that u(t,x) = Y t ' x is the unique (in Ci 3 ) 

stochastic viscosity solution of SPDE (4.9). This extends the Feynraan-Kac formula to SPDEs driven 

by a fractional Brownian motion. 

We conclude the main theorems of Section 4 with the following relation diagram, which shows the 
mutual relationship between fractional backward SDEs and SPDEs: 



u(t, x) is the viscosity solution of (4.8) 

X 

(Y, Z) is the solution of BSDE (4.5) 



GT 



GT 



u(t,x) is the viscosity solution of (4.9) 

X 

(Y, Z) is the solution of BDSDE (4.4) 



where 'GT' stands for 'Girsanov transformation'. 

Finally, in order to illustrate how our method works, we give the example of a linear fractional 
backward doubly stochastic differential equation. 

Example 4.20. We let d = 1 and f(s 1 x 1 y,z) = fjx + fly + f^z, where the coefficients f},fg 
and fg are bounded and deterministic functions. The associated fractional backward doubly stochastic 
differential equation is linear and writes: 



Y*> x - $(X*' X ) + / (flXl' x + f r Y^ x + f r ZY) dr - / Zl' x i dW r 
Jo Ja 



lrY*> x AB, 



(4.29) 
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After Girsanov transformation, it becomes 

Y^ = $(X*' X ) + J ^X* r ' x e- 1 (T r ) + f 2 r Y l r - x + flZ 1 ^ dr - f Zp x 1 dW r . 
and has the following solution: 

Y s l ' x = E Q J' (jlX^s-^Tr) cxp y S tiduj^j dr + exp j jT / r 2 dr| $(^' x ) 

where Eq is the expectation with respect to Q — cxp i J Q f^dW r — 5 Jo (fr) 2( ^ r \ P- According to 
orem 3.9, the solution of (4.29) is then 



^vj-f 



The- 



Yf'* - E Q 



J (f^X^e; l (T r A s ) cxp U f^duX) dr + e s cxp U f r dr\ $(X, 



O' / 



^ V Tf 
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